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Abstract: 0 4 Klein Gordon Equation has been solved numerically by using two methods: finite difference 
method (FDM) and Adomain decomposition method (ADM) and we discover that the ADM is much more 
accurate than FDM in this kind of models as shown in the examplef 1,2). 
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I. INTRODUCTION 

A finite difference method proceeds by replacing the derivatives in the differential equations by finite 
difference approximations. We assume the function whose derivatives are to be approximated by using Taylor's 
expansion. This gives a large algebraic system of equations to be solved in place of the differential equation, 
something that is easily solved on a computer [4]. 

The main advantage of the method is that it can provide analytical or an approximated solution to a 
rather wide class of nonlinear (and stochastic) equations without linearization, perturbation, closure 
approximation, or discretization methods. Unlike the common methods which are only applicable to systems 
with weak nonlinearity and small perturbation and may change the physics of the problem due to simplification, 
ADM gives the approximated solution of the problem without any simplification. Another important advantage 
is that the method is capable of greatly reducing the size of computation work while still maintaining high 
accuracy of the numerical solution [1,2]. 

1.1. MATHEMATICAL MODEL: 

The Klein-Gordon equation, 

g = c 2 Au-/(u) (1) 

was named after the physicists Oskar Klein and Walter Gordon, who in 1926 proposed that it describes 
relativistic electrons. Some other authors make the similar claims in the same year. The equation was first 
considered as a quantum wave equation by Schrodinger in his search for an equation describing de Broglie 
waves. The equation is found in his notebooks from late 1925, before he made the discovery of the equation that 
now bears his name. He rejected it because he couldn't make it fit data (the equation doesn't take in account the 
spin of the electron), the way he found his equation was by making simplifications in the Klein Gordon 
equation. In 1927, soon after the Schrodinger equation was introduced, Vladimir Fock wrote an article about its 
generalization for the case of magnetic fields, where forces were dependent on velocity, and independently 
derived this equation. Both Klein and Fock used Kaluza and Klein's method [6]. 

When/(u) = sin (u), then equation (1) becomes sine-Gordon equation, which is found by Zabusky 
and kruskal in 1965. 

Fiore et al. (2005) gave arguments for the existence of exact travelling wave solutions of a perturbed sine 
Gordon equation on the real line or on the circle and classified them [7 ]. 

When f(u) = mu — £U 3 then equation (1) is called 0 4 -nonlinear klein Gordon equation (0 4 
equation)[8] 

— = c 2 — - mu + eu 2 (2) 
or [33] 

d 2 u d 2 u 2 , 0 i 

With initial and boundary conditions [10] 
u(x, 0) = f(x) and = 0 ,t > 0 ,0<x<2n 

u(0, t) = u(2n, t) = 0 

Equation (3) arises in quantum field theory with m denoting mass and e is 
coupling constant [11]. 

The <p 4 equation has become an important subject because of its numerous applications in condensed 
matter physics. It describes, for example, structural phase transitions in ferroelectric and ferromagnetic 
materials, topological excitations in quasi one dimensional system like biological macromolecules and hydrogen 
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chains, or polymers, etc. Its simplest localized solutions so-called "kinks" which are related to the motion of the 
aforementioned topological excitations, e.g., domain walls in second order phase transitions, or polymerization 
mismatches. A more realistic modeling of physical situation in condensed matter physics often requires the 
inclusion of perturbations of different types like thermal noise and time or spatial dependent potential 
fluctuations [ 12 ]. 

II. MATERIALS AND METHODS 

II.1 Finite Difference Method 

Our aim is to approximate solutions to differential equations (i.e., to find a function or some discrete 
approximation to this function) which satisfies a given relationship between various of its derivatives on some 
given region of space and/or time, along with some boundary conditions along the edges of this domain. Now, 
the formulation of this technique gives: 

, = .N)-W 

v ' h 

u"(x) = — (u(x + h) — 2u{x) + u(x — K)) 
Let the coordinate (x, t) of the grid points be 

x = ph, t — qk 

Where p, q are integers. 

Denote the values of u at these mesh points by 

u(ph,qk) = u pq 

We start by partitioning the rectangle R = {(x, t):0<x< 2n, 0 < t < b} into a grid consisting of n — 1 by 
m — 1 rectangle with sides Ax — h andAt = k, as shown in Figure (1) [4]: 
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Figure (1): The mesh of finite difference method 

We shall use a difference-equation method to calculate approximation 

{u pq : p = 1,2,3, ...,n} ,q — 1,2,3,4, ...,m 

Now, equation (3) be as follows: 

u V,q+l- 2u V,q+ u V,q-l _ u v+l,q- 2u v,q+ u p-l,q , 3 

fc 2 h 2 " tu p,q ' K u p,q 

u (p,q+i) ~ 2u p q + Up q^ = — 9 ) — 2u p q + u p _ l q ) — mk 2 u p q + £ji 2 U p cl 

k 2 

setting r — — 
h 

u p,q + l ~ ^Up q + U p q_ 1 = r(u( p + 1 ^ — 2u p q + U p _ l q ^ — mk 2 U p q + £ji 2 U p q 
U p,q + 1 ~ r (_ U (p + l,q) + U p-l,q) ~ ^ rU p,q + 2lL p q — U p q_ 1 — mk 2 U p q + eJc 2 U p q 
U p,q + 1 = r ( U (p + l,q) + U p-l,q) + (2 - 2f - mk 2 + ^U^Upq - Up,^ (4) 

Equation (4) represents the explicit difference approximation for 0 4 equation, which is used to find the 
row (q + 1) across the grid where the approximation in both rows (q), (q — 1) are known. The four known 
values on the right hand side U( p U(p-i ,q> u (p+i,q) an( J u (p,q) are use d t0 create the approximation U( p „ +: n . 
In order to use formula (4) to compute the third row, the first two starting rows of values corresponding to q = 1 
and q — 2 must be supplied. The boundary function f{x pi t x ) = f(x p ), p = 1,2,3, n is used to help 



www.ijeijoumal.com 



Page I 38 



Adomain Decomposition Method for 0 4 Klein Gordon Equation 



produce the first row. Since the second row cannot be determined in a usually way, we may use the Taylor's 
formula u{x, t) at x = x p to compute it. 

Up,2 = Up,i + r -(u p+lil +u p _ lil ) - ru Pil - Y mu p,i + \ £U l,i ^ 

II.l Adomain Decomposition Method 

We will present a review of the standard Adomin decomposition method [1,2,3]. To achieve our goal 
we consider the general form of differential equation: 
Lu + Ru + Nu — g (x) 

Where L is the highest order derivative which is assumed to be easily invertible R is the remainder of 
the linear operator, Nu represents the nonlinear terms, and g is the source term. 
The Adomin's technique consists of approximating the solution an infinite series: 

U(x,t) =£n=o"n (X,t) 

And decomposing the nonlinear operator N as 

N(u)=Y,% =0 A n 
Where A n is Adomian's polynomial of u 0 , u v u 2 ,u 3 ,u n 

A n = ^[N(lT=o^ t )]^o,n>0 

and we consider the nonlinear homogenous Klein-Gordon 0 4 equation: 

u tt - u xx + mu ~ £V - 3 

Applying the decomposition method (see [1], [2], [11]) 

L tt u = L xx u + mu - eu 3 

Where L tt = ^ and L xx = ^ 

Assuming the inverse of operator L tt exist, it can be taken as 

rt = /o/oC) dtdt - 
Therefore, applying on both sides with LJt yields 

u(x, t) = u(x, 0) + tu\x, 0) + Ltt(L xx u) + mhllu — eL^u 3 
By using initial condition we get 

u(x, t) = u 0 (x) + tu 0 (x) + L~ll(L xx u) + mLT t tu - eL^u 3 

So, we get: 

Yn= 0 u(t, x) = u 0 (X) + tu 0 {x) + L- t £{L xx Y£ =0 u) + mL- t £ J% =0 u - £L- t £ Xn=0 K 
Then we have: 
u 0 = u 0 (x), 
and 

"fc+i = LTt<iLxxU k ) + mLttW - £L~ t lA k 

Where k > 0. 
For k = 0 , 

A 0 = ~[F(Zl 0 A°u 0 -)] A=0 = F(u 0 ) = u 3 

Ao — u o 
Then we get , 

u ± = L^{L xx u 0 ) + mL^u 0 - eL^u 3 , = L^iL^u^t + ml^u^t - sL^u^t 
= ((LxxU-o) + mu 0 - eu 3 ,) - = {L xx u 0 ) - + mu 0 - - eu% -, 

Let Dl = ((L xx u 0 ~) + mu 0 - eu 3 ,) 
t 2 

Then u-, — Dl — 

1 2 

For k — 1 

A ± = jij^[F(Xl =0 A i u i )] x=0 =^-[F(u 0 +/l 1 u 1 )] A= o=^-[(u 0 +/L 1 u 1 ) 3 ] A= o[3(u 0 +/l 1 u 1 ) 2 u 1 ] /l= o=3ugu 1 
we get: 

u 2 = ^(inUi) + mLttV - £Ltt 3uo"i = Lt 1 (L xx u 1 )t + ml^u^t - 3sLt 1 u%u 1 t 

= 7 + mu ± — - 3euIu ± - = {l xx u v + mu ± - 3eu^u 1 ) — 

Let, D2= + - 3^0%) 

Then u 2 = D2- 

Z 2 

For fe = 2 

A 2 = ^[FGLo^i)]^ =^K u o + A 2 u 2 ) 3 ] A=0 =^[(u 0 +AV +A 2 u 2 ) 2 ( Wl + 2Au 2 )] x=0 
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=^[( u o + A 2 u 2 ) 2 (2u 2 ) + 2(u 0 +A 1 u 1 + A 2 u 2 )(u 1 + 2Au 2 )(u 1 + 2Au 2 )] A=0 

=3u 2 u 2 + 3u 0 u 2 

We get: 

u 3 = Ltt(L xx u 2 ) + mL~llu 2 — £L7 £ 1 (3u 2 u 2 + 3u 0 u 2 ) = Ll 1 {L xx u 2 )t + mLt 1 u 2 t — £L7 1 (3u 2 u 2 + 3u 0 u 2 )t , 



= L xx (u 2 ) j + mu 2 j - e(3u 2 u 2 + 3u 0 u 2 ) l — D3 = L xx (u 2 ) + mu 2 - e(3u 2 u 2 + 3u 0 u 2 ), 

= US—, 

6 2 

and so on for k — 3, 4, ... 

III. APPLICATION (NUMERICAL EXAMPLES) 

We solve the following examples numerically to illustrate the efficiency of the presented methods, suppose we 
have the system 



Example 1: [6] 

d 2 u d 2 u 



— 7 = —a —r — Bu — yu 2 

dt 2 dx 2 r ' 



,t > 0 

With the initial conditions 

u(x, 0) = Btan(Kx),u t = BcK sec 2 (Kx) 

We take a = -2.5, B = 1, y = 1.5. 

Where |/? 



-1 <x < 1 




Fig. (2) ADM solution to 0<t<0.01 
and -1<x<1 with dx=0. 105263 157894737 




Fig. (3) EFDM solution to 0<t<0.01 and -1<x<1 
with dx=0. 105263 157894737 




Fig.(4)Exact solution to 0<t<0.01 and -1<x<1 with dx=0. 105263 157894737 
Table (1) comparison between ADM, EFDM and Exact solutions at t=2 with dx=0.020408 163265306 



EXACT 


ADM 


EFDM 


-0.414411214133078 


-0.414411214101644 


-0.414411214133078 


-0.364670241057009 


-0.364670241031911 


-0.364669012962060 


-0.317087953610942 


-0.317087953590899 


-0.317086928637545 


-0.271306324156148 


-0.271306324140212 


-0.271305477633181 


-0.22701112796834 


-0.227011127955811 


-0.227010440389676 


-0.183922496540239 


-0.183922496530586 


-0.183921952741587 


-0.141787250668732 


-0.141787250661578 


-0.141786839088029 


-0.100372541039790 


-0.100372541034864 


-0.100372253179315 


-0.059460438539449 


-0.059460438536567 


-0.059460268585494 


-0.018843195109709 


-0.018843195108769 


-0.018843139679920 
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0.021681051104011 


0.021681051103044 


0.021680993104744 


0.062312318757989 


0.062312318755080 


0.062312146183729 


0.103252748218675 


0.103252748213720 


0.103252457634093 


0.144710656400331 


0.144710656393147 


0.144710241934291 


0.186904856153030 


0.186904856143345 


0.186904309243295 


0.230069434484844 


0.230069434472276 


0.230068743494559 


0.274459218894519 


0.274459218878543 


0.274458368570347 


0.320356215567061 


0.320356215546973 


0.320355186294337 


0.368077383902961 


0.368077383877811 


0.368076150876117 


0.417984229397421 


0.417984229365926 


0.417984229397421 



Example 2: [see (6)] 

d 2 u d 2 u o _ 

— - = — - — u — u* ,t > 0 

dt 2 dx 2 

With the initial conditions 

u(x, 0) = A[l + cos (^)] , u t = 0 , 0 < x < L 
The boundary conditions are given by 
u x (0, t) = 0, u x (L, t) = 0 , where L=l and A=1.5 





Fig. (5) ) ADM solution to 0<t<0.0001 
and 0<x<l with dx=0.020408163265306 



Fig. (6) EFDM solution to 0<t<0.0001 
and 0<x<l with dx=0.020408163265306 



Table (2) compari son between ADM, EFDM solution at t=2 with dx=0. 020408 163265306 



ADM 


EFDM 


2.987685020734869 


2.987685020550971 


2.950942294558544 


2.950942294379764 


2.806978056185084 


2.806978056025879 


2.577524025146591 


2.577524025017099 


2.435234702788101 


2.435234702676172 


2.107175014683591 


2.107175014609870 


1.739399842550069 


1.739399842515833 


1.548077366357483 


1.548077366342586 


1.355965461138478 


1.355965461142287 


1.166218599065529 


1.166218599087200 


0.981952418368039 


0.981952418406554 


0.806192564638747 


0.806192564692939 


0.641825009816746 


0.641825009885315 


0.491548664608025 


0.491548664689548 


0.357831062446298 


0.357831062539230 


0.242867842662239 


0.242867842764914 


0.148546698146371 


0.148546698257009 


0.076416379483997 


0.076416379600713 
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O.Oz/oo 1 2645 1 3402 


0.0z7oolz64o34zz3 


A A AT AO 1 A1 AOT/I /I AC 

0.00308 191 0874495 


A A AT AO 1 A 1 A A AT} O/C 

0.003081910997386 


A A AO AO 1 A1 AOT/I ,(AC 

0.00308 191 0874495 


A A AO AO 1 A 1 AAmO O tl 

0.003081910997386 


U.Uz/oo 1 Zo45 1 340Z 


0.0z7oolzo4o34zz3 


U.U / 0410 J /y4oj77 / 


U.U/0410j/!70UU/ 1j 


U.14oJ400yo 140.5 / 1 


U. 14oj400VoZj /UUV 


U.Z4Z00 / 04Z00ZZ 


M O/lOQATC/lO^A/lQl /I 
U.Z4ZoO /o4Z /04V 14 


V. J J 1 oj 1U0Z440ZV0 


U.jj /ojlVOZjjyZZy 


U.4V 1 j4o0040UoUZ4 


U.4V 1 j4o0040oVj4 / 


a / ( i C AAAO \ tl'l A d 

0.64 1 8250098 1 6745 


A £/| 1 O^C AAA OOC11 A 

0.641825009885314 


0.806192564638747 


0.806192564692939 


0.981952418368038 


0.981952418406553 


i . luuziojyyuojj/ / 


i . loo/i ojyyuo / iyo 


1.355965461138476 


1.355965461142286 


1.548077366357482 


1.548077366342585 


1.739399842550069 


1.739399842515833 


1.926791379946547 


1.926791379892598 


2.107175014683589 


2.107175014609868 



IV. CONCLUSION 

We saw that Adomain decomposition method is much more accurate than finite difference method for solving 
□ 4 Klein Gordon Equation and this kind of models as shown in figures (2-4) and table (1) for example (1) and 
figures (5-6) and table (2). 
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